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Abstract 

The influence of the spin-orbit couphng on the magnetic structure of deposited transition metal 
nanostructure systems has been studied by fully relativistic electronic structure calculations. The 
interplay of exchange coupling and magnetic anisotropy was monitored by studying the correspond- 
ing magnetic torque calculated within ab-initio and model approaches. It is found that a spin-orbit 
induced Dzyaloshinski-Moriya interaction can stabilise a non-collinear spin structure even if there 
is a pronounced isotropic ferromagnetic exchange interaction between the magnetic atoms. 
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As the complexity and technological applicability of nanostructured magnetic materials 
grows it has become important to develop a reliable quantitative theoretical framework in 
which to understand them. In principle this is available from relativistic density functional 
theory (RDFT) []]|. The magnetic behaviour of complex magnetic systems is described 
effectively in terms of 'local moments', even in metallic systems, which fluctuate as the 
temperature is increased from OK. Within the effective single electron picture of DFT, the 
moments arise from local magnetic fields associated with the different atomic sites which 
affect the electronic motions and are self-consistently maintained by them. The spin-orbit 
coupling (SOC) effects on the electronic structure determine the magnetic anisotropy linking 
the magnetic and spatial structure of a material. As demonstrated by Staunton et al. 2| 
and others, magnetic torque calculations enable magnetic anisotropy to be studied reliably 
and magnetic structures to be determined. We have recently developed a method to study 
ab-initio the magnetic structures of complex nanostructures using this approach sl. It is 
based on the framework of relativistic density functional theory using the local spin density 
approximation (LSDA) for exchange and correlation effects j4|. The electronic structure was 
determined in a fully relativistic way from the basis of the Dirac equation for spin-polarized 
potentials which is solved using the Korringa-Kohn-Rostoker (KKR) multiple scattering 
formalism j^. 

In this letter we set out to understand the relativistic effects on the magnetic structure of 
nanoclusters in simpler terms of a model 'moment' Hamiltonian. In this context the map- 
ping of energetic properties obtained from first principles calculations of complex magnetic 
systems onto a Heisenberg Hamiltonian has proved over recent years to be a very robust 
and successful scheme g], Q]. Moreover for many systems the results concur with those from 
an ab-initio 'disordered local moment' theory in which no prior mapping to a Heisenberg 
system is assumed . An extensively used approach to calculate the isotropic exchange 
interaction parameter Jij for two magnetic moments on sites i and j for use in the classical 
Heisenberg Hamiltonian was worked out by Lichtenstein et al. Gj using perturbation the- 
ory and the so-called Lloyd formula. A corresponding fully relativistic approach was later 
introduced by Udvardi et al. [7] that produces an exchange interaction tensor J_._. for use 
in a prescribed extended Heisenberg Hamiltonian. This scheme generates in particular a 



Dzyaloshinski-Moriya (DM)-type interaction [lO, lUj, which may explain many interesting 
phenomena including the magnetic ground state configuration of nanostructures [12I, Q as 
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well as maEfnetic thin films , [l5| . Here we explore what form the effective moment Hamil- 
tonian should have by a detailed study of its building blocks. With an investigation of Fe2, 
C02 and Ni2 dimers deposited on Pt(lll) as most simple cluster examples we identify strong 
DM-type interactions as well as an additional substrate-generated effect in Ni2 (for calcula- 
tional details see ^]). The conclusions of these results can then be applied straightforwardly 
to larger clusters and nanostructures. 

We start by considering the magnetic torque vector Tj^*^^^ acting on an atomic magnetic 
moment on a site i {i = 1 or 2) and aligned along direction Cj. The torque vector is defined 
in terms of the change in energy E{{ek}) = E{ei,e2) of the system when changing the 
orientation of the magnetic moment, on site i, Tj^^"^ = —dE{{ek]) / dci. The component 
= -{dE{{ek})/dei) ■ {u x e^) of ff with respect to the axis u can be determined from 
first-principles using an expression derived by Staunton et al. Following on from this 
the derivative c^EI3,.9e, deser.biug the change ■„ energy upon cha„g,„g the onentation 
of two magnetic moments on sites i and j can also be obtained [6|, 17|. By making use of 
the rigid spin approximation (RSA) [l6| this approach leads to a fitting of the magnetic 
energy landscape E{{ek}) at low temperatures obtained from first principles calculations 
onto a Heisenberg model Hamiltonian. With SOC included an anisotropy in the exchange 
interaction may occur. We then use these quantities in an extended classical Heisenberg 
Hamiltonian for the 'spins', {cfc} of the following form, (e.g. 

,0): 



-\ Di,-[e,xe,] + Y,Ki{e,) ■ (1) 

Here the exchange interaction tensor J_.^ has been split into its conventional isotropic part 
Jij, its traceless symmetric part J_^. and its anti-symmetric part J^. We assume that the 
latter one is represented in terms of the DM vector Dij with D]j = e°^^ '-^ ^ (^^°'f^'y jg 
the Levi-Civita symbol). Finally, the anisotropy constants Ki{ei) account for the so-called 
on-site anisotropy energy associated with each individual moment oriented along Cj. From 
our ab-initio calculation of the torques T^*'^"'' we can test whether such a model Hamiltonian 
is justified and also find the values of the Jij, Ki, J_^., and Dij parameters. 

For the model Eq. ([1]) the rate of change in energy when a moment on site i is rotated 
about an axis u can be partitioned into contributions from different terms of Heisenberg 



Hamiltonian T.J^ = T}^ + T^,u + T^^i^ + T^^u- The contribution to the torque from the DM 
couphng is given by: 

The last term derives from the single site anisotropy term, i.e. = ^^q-'^^^ • {u x Cj). 
By focussing on nanoclusters of collinear 'spin' arrangements, the magnetically anisotropic 
terms can be determined. Note that the last term in Eq. ([2]) does not contribute to the 
torque in the case of a collinear magnetic structure. Also the sum of the DM contributions 
from all sites in a nanocluster to the total torque vanishes in this case. On the other 
hand, the anisotropy of the exchange interaction, represented by the symmetric tensor J^. 
gives rise to a finite torque hence contributing to the total magnetic anisotropy energy of 
the nanocluster. 

We now see how well our ab-initio calculations of deposited Fe, Co and Ni dimers, using 
our first-principles multiple scattering formalism, fit such model torques. Firstly we find 
that the dimers, when constrained to be magnetically collinear exhibit a pronounced out- 
of-p,a,. .ag,.t,c a„,sot.opy R, The co,.,„ea. magnetic moments Ci and 62 have been 
assumed to be orientated at an angle to the surface normal (z-axis) and are expressed in 
terms of polar angles 6 and (p, as ei = 62 = (sin 6' cos 0, sin 6' sin 0, cos ^). Fig. [1] shows the 
atomic configuration together the projection of the moments onto the surface {xy-p\ane) . 
In parallel to the fixed frame of reference (x, y, z) we use a second one (x', y' , z') rotated by 
with respect to the fixed one with z = z' (see Fig. [1]). In all calculations the torque T^^^ 
is taken around the y' axis, i.e. u = y' = {— sin0, cos0, 0), for 9 fixed to 7r/4 as a function 
of the azimuthal angle 0. From the model Eq. ([T]) an expansion of the anisotropy energy 
term K{ei) in terms of Z = 2 spherical harmonics, gives the single site contribution to the 
torque in present configuration to be = —2[K2^i + -ft'2,2 cos(20) + i^'j 2 sin(20)] whereas 
T^^^ deduced from the model has a cos(0) (sin(0)) variation when D12 lies in the yz (xz) 
plane. 

In Fig.[2]the ab-initio results for the torque Tf" are compared with those deduced from the 
model Eq. ([T]). As one notes, the torques for the two atoms are different but are related with 
respect to their 0-dependence according to the symmetry of the system. On comparing 
the torque of the two Co atoms of Co2/Pt(lll) obtained directly from the electronic structure 
calculations with those from the model (see Fig. [2]a-c) we find the symmetric part of the 
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FIG. 1: Magnetic configuration of the transition metal dimers deposited on a Pt(lll) substrate. 
Tlie large (small) spheres present Pt surface (subsurface) atoms. The medium size spheres repre- 
sent the dimer atoms with the projection of their magnetic moments onto the surface (j;y-plane) 
represented by arrows. The projection of the DM-vector D12 onto the surface is represented by a 
short arrow. 

exchange interaction tensor to be negligible. The contribution Tf'^ due to the DM 
interaction Eq. ([2]) is shown in Fig. [2f). Clearly, and vary with coscp and are 

opposite in sign. In accordance with the symmetry, the DM vector D12 lies in the yz 
plane (Fig. [1]). For the contribution T^^ due to the on-site anisotropy one finds K'22 to 
be very small and the dominating terms i^2,i and K2^2 to be practically the same for both 
atomic sites leading to ^ T^. 

The contribution to T^^ connected with K2^i does not depend on 0, while that connected 
with K2^2 varies with cos(20) (Fig. [2]c). As one can see, the torque T™°'^ derived from the 
model Hamiltonian reproduces the results T^^^ calculated directly rather well. The remaining 
deviation is primarily due to the limitations of the model Hamiltonian with respect to the 
dependency of magnetic energy on the magnetic moment orientations E{{ek})- From the 
decomposition of the Tj's via the model Hamiltonian it becomes clear that its ^-dependence 
is dominated by the DM-contribution while the /r2,2-contribution gives rise to a minor 
additional modulation. Owing to the large negative values of Ti{K2^i) and T2{K2,i) an out- 
of-plane anisotropy results for the total system. This can be seen also from Fig. [2]e where 
the total torque T^^l = T^" + T^"' for C02 on Pt(lll) is shown together with the individual 
contributions Tf^"^ and the corresponding DM terms T^^. Once more one sees that the 
0-dependence of the individual torques is determined by while that of the total torque 
T^^l is set by the 7^^2,2 on-site anisotropy terms. This also holds for Fe2 (Fig. [2H) for which 
the DM terms are even more dominant, i.e. the 0-dependence of the individual torques are 
nearly exclusively due to TP^. 
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FIG. 2: Magnetic torque Tj for dimers on Pt(lll) for = 7r/4 as a function of the azimuthal angle 
(p. Left: Torque for Co2/Pt(lll): a) results Tf" of direct calculations (up and down triangles) 
compared with results T™°*^ obtained on the basis of the Heisenberg-type model Hamiltonian 
(Eq. dH)); 6) T^^ contribution to Tf^""^ due to the DM interaction (Eq. ©); c) contribution Tf^ 
to Tf^°'^ due to the on-site anisotropy. Right: Magnetic torque for the two atoms of Fe2 (d), 
Co2 (e) and Ni2 (/). The results calculated directly are presented by up (T{*"') and down (T^"') 
triangles and their sum by thick solid line {T^l = T{^"' + T^''^). The thin solid lines give the DM 
contributions according to Eq. ([2|). 

The situation is noticeably different for Ni2 on Pt(lll) (Fig. [2]f) for which the DM terms 
give only minor contributions to the individual torques Ti and T2. In contrast to Fe2 and C02 
the difference between the torques of different atoms cannot be attributed to DM coupling. 
The period of oscillations of the torques as functions of (p are different compared to those 
of the torques created by DM coupling varying as sin(20). Clearly an effect which is not 
contained in the model Hamiltonian Eq. ([1])) is evident and must derive from interactions 
of the Ni moments with those induced in the Pt substrate. As with a DM-type interaction 
it produces torques such that Ti = — T2 but its sin(20) variation would arise if a term of 
the form {A ■ ei){B ■ ei) — {A ■ 62) ■ {B ■ 62) were added to Eq. ([1]), where A lies along x 
and B along y. These Ni2 dimer results indicate that an effective Heisenberg model must 
be used with caution for systems where the magnetic structure of a nanocluster is strongly 
influenced by the spin polarisability of the substrate. 

For the C02 and Fe2 dimers, however, the model works very well. For our chosen geometry 
6 = it/A and = and using the symmetry properties of the elements of the exchange tensor, 
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TABLE I: Components of the DM vector Dij, the isotropic exchange constant Jij (in meV) 
and the tilt angle a (in degrees) (see text) for the dimers Fe2 and C02 on Pt(lll). The data 
labelled direct have been obtained from the direct first principles calculations of the torque. The 
data below have been obtained by mapping first principles results onto the model Hamiltonian, 
Eq. dH). 





direct 


model 








^12 




^12 


J12 


a 


Fea 


6.04 


0.00 


6.07 


-3.34 


138.0 


2.52 


C02 


3.69 


0.00 


3.89 


-3.84 


108.0 


2.07 


Ni2 


-0.02 


0.00 


-0.04 


-0.24 


30.4 


0.07 



one finds for Co2/Pt(lll) with ^ the total torque T = Ti + T2 = -( Jfs"^ - ^1^2"'") + 
2rf . Finally, the magnetic anisotropy energy (MAE) of the dimer being the difference in 
energy when the magnetic moments are both oriented along and e^, is given by the integral 
— j^^{Ti{e) +T2{e))de js]. Obviously this has no contribution from the DM interaction. For 
Co2/Pt(lll) we find the exchange parameters Jj^^, Jf| and Jf| to be nearly the same 
implying that the total MAE of the dimer is nearly exclusively due to the on-site anisotropy. 
The values i^2,i = 1-5 meV and i^2,2 = 0.39 meV for Co2/Pt(lll) lead as mentioned above 
to a pronounced out-of plane anisotropy, i.e. in the ground state the total magnetisation 
points along the surface normal. 

Taking the difference between the individual torques one arrives at the relation Df^ = 
^J^^ allowing be deduced directly from the ab-initio torque Tf^^ calculations. Table 

[T] shows the corresponding results for dimers Fe2 and C02 on Ft (111) in comparison with 
data derived from a mapping to our model Hamiltonian, Eq. ([T]). 

The closeness of the results justifies once more the use of the model Hamiltonian for Fe2 
and C02. In addition one notes that 0^2 has an appreciable value compared to the isotropic 
exchange constant J12. Fixing the azimuthal angle to be 7r/2 and performing similar steps 
one finds Df2 to be 0. This is also in line with the Cg symmetry of the investigated dimer 
systems (see Fig. [1]). DI2, on the other hand, may take a non-zero value and is found to be 
comparable to D\2 (see Table H]). 
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Thus, the above analysis shows that the torques Ti and T2 may differ even if the total 
torque is zero, i.e. if the moments are aligned coUinearly along the easy axis (surface normal). 
The difference between these torques is caused exclusively by the terms leading to a 
rotation around the y-axis. Minimizing the magnetic energy E{{ei}) of the 2 atom clusters 
leads to an outward tilting of the moments by an angle a given by a = atan (-D12/ "^iz) 
with J12 being the isotropic exchange coupling constant. The corresponding results given in 
Table [T] show that the DM interaction causes the deposited Fe2 and C02 dimers to have an 
appreciable deviation from collinear configurations in spite of the pronounced ferromagnetic 
exchange coupling given by J12. This effect of SOC is completely in line with the findings 



of Sandratskii and Kiibler 



19| for bulk systems. 



The substrate clearly plays a crucial role in the DM interaction. Firstly, hybridisation 
with the substrate breaks the inversion symmetry for the dimer leading to a non-zero DM 
vector. This symmetry effect is also confirmed by test calculations on free dimers and on 
dimers embedded in a bulk Pt matrix. Secondly, the hybridisation with the substrate also 
allows the SOC effects of the substrate to be transferred to the magnetic 3d transition metal 
dimer. We have confirmed this by further calculations in which the SOC of the substrate 
and the dimer atoms were manipulated separately. Enhancing the SOC for a Co dimer leads 
primarily to an increase of the on-site anisotropy -ft'2,2- However, enhancing the SOC for the 
Pt substrate leads to a strong increase in the anisotropy 7^2,1 as well as to a larger difference 
in the individual torques on the two Co atoms, reflecting an increase of the DM interaction 
(Fig. [2]). This behaviour is in line with Levi's model of the indirect DM interaction between 
two spin moments 201], which is mediated by nearby atoms. As a consequence, the magnitude 
of the DM interaction is essentially determined by the SOC strength of the neighbouring 
atoms. 

The effect of anisotropic exchange is even more spectacular in magnetic alloy nanoclus- 
ters, where magnetic atoms are separated by non-magnetic atoms with large SOC. This is 
demonstrated by our calculations for a FePt (2 x 1) 2D alloy cluster deposited on the Pt(lll) 
surface pertinent to an experimental investigation by Honolka et al. 21[. Ab-initio torque 
calculations show the clear trend that the magnetic structure of this cluster is non-collinear. 
Here we present the results of Monte Carlo simulation based on an effective Heisenberg 
model with the isotropic and DM exchange coupling parameters calculated in accordance 
to Ref. ?!. Fig. [3] shows the non-collinear magnetic structure obtained for T = OK. The 
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FIG. 3: Magnetic structure of a 91-atom FePt (2 x 1) 2D alloy cluster deposited on the Pt(lll) 
surface determined via MC simulation for T = OK. The arrows denote the orientation of the 
magnetic moments associated with the iron atoms. 

non-collinearity between the Fe chains is essentially caused by the nearest neighbours Fe-Fe 
interchain DM interaction {\D\ = 4.6 meV) being of similar magnitude when compared to 
the isotropic exchange interaction (J = 8.8 meV). Within a Fe chain, however, the DM 
interaction is more than one order of magnitude smaller when compared to the isotropic 
exchange leading only to a slight screwing of the Fe magnetic moments along the chain. 
Thus, we see that the FePt cluster's non-collinear structure is created by an enhanced DM 
coupling between Fe moments mediated by Pt atoms having large SOC, on the one hand 
side, while the isotropic exchange between these atoms separated by non-magnetic Pt is 
small enough to make both these couplings comparable. 

In summary our investigation of deposited transition metal dimers as simple but realistic 
model systems has demonstrated that ab-initio magnetic torque calculations enable the im- 
pact of SOC on the magnetic interactions within nanostructures to be monitored in a very 
detailed way revealing subtle anisotropic effects. The analysis of directly calculated elec- 
tronic structure quantities within a framework based on a Heisenberg model Hamiltonian 
gives further insight and identifies the role of the various contributions and also the limita- 
tions of such models. For the Fe and Co dimer systems studied here the DM interaction was 
found to be pronounced owing primarily to the SOC of the substrate. Moreover it leads to 
non-collinear magnetic configurations of the dimers in spite of the pronounced ferromagnetic 
coupling and out-of-plane anisotropy. These magnetically anisotropic interactions also have 
a profound effect in larger clusters and we have demonstrated this explicitly with a study 
of a deposited FePt cluster. In particular we infer that the magnetic structure around the 
edges of magnetic nanoparticles is likely to be significantly affected. 
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